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HIGHLIGHTS 


•  A  simple  but  accurate  battery  model  based  on  equivalent-circuit  is  required. 

•  Multi-objective  genetic  algorithm  is  utilized  for  extracting  performance  parameters. 

•  Battery  performance  is  accurately  predicted  once  parameters  are  optimally  extracted. 

•  The  parameter-extracting  code  is  widely  applicable  for  various  types  of  batteries. 

•  It  can  serve  as  a  robust  and  reliable  tool  for  extracting  the  battery  parameters. 
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A  simple  but  reasonably  accurate  battery  model  is  required  for  simulating  the  performance  of  electrical 
systems  that  employ  a  battery  for  example  an  electric  vehicle,  as  well  as  for  investigating  their  potential 
as  an  energy  storage  device.  In  this  paper,  a  relatively  simple  equivalent  circuit  based  model  is  employed 
for  modeling  the  performance  of  a  battery.  A  computer  code  utilizing  a  multi-objective  genetic  algorithm 
is  developed  for  the  purpose  of  extracting  the  battery  performance  parameters.  The  code  is  applied  to 
several  existing  industrial  batteries  as  well  as  to  two  recently  proposed  high  performance  batteries 
which  are  currently  in  early  research  and  development  stage.  The  results  demonstrate  that  with  the 
optimally  extracted  performance  parameters,  the  equivalent  circuit  based  battery  model  can  accurately 
predict  the  performance  of  various  batteries  of  different  sizes,  capacities,  and  materials.  Several  test  cases 
demonstrate  that  the  multi-objective  genetic  algorithm  can  serve  as  a  robust  and  reliable  tool  for 
extracting  the  battery  performance  parameters. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  recent  years,  there  has  been  tremendous  emphasis  as  well  as 
effort  devoted  towards  the  development  of  high  performance 
batteries  mainly  driven  by  their  need  as  a  power  source  for  an 
electric  or  hybrid  electric  vehicle  [1]  and  as  an  energy  storage  de¬ 
vice  for  intermittent  energy  supply  sources  such  as  solar  and  wind 
[2],  Rechargeable  batteries  are  particularly  suitable  for  these  ap¬ 
plications  because  of  their  potential  for  high  round-trip  efficiency, 
long  operating  life  and  scalability  [3],  In  addition  to  these  technical 
requirements,  the  battery  should  be  of  low  enough  cost  to  be  cost 
effective  as  a  power  source  for  an  electric  vehicle  or  as  an  energy 
storage  device. 
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To  study  the  behavior  of  electric  vehicles,  manufacturers  model 
the  electric  drive  train.  In  order  to  do  so,  the  energy  source  of  these 
drive  trains,  the  battery,  needs  to  be  accurately  modeled  as  well. 
There  are  essentially  four  types  of  battery  models  that  are  used  to 
model  its  performance:  these  can  be  categorized  as  the  experi¬ 
mental  models,  the  electrochemical  models,  the  mathematical 
models  and  the  electric  circuit  based  models.  The  experimental  and 
electrochemical  models  are  not  well  suited  to  model  the  cell  dy¬ 
namics:  hence  they  cannot  be  used  reliably  to  determine  the  State- 
of-Charge  (SOC)  of  the  battery  [4],  Furthermore,  the  electro¬ 
chemical  models  are  computationally  expensive  and  require 
extensive  experimentation  for  determining  the  parameters  of  the 
model  [1],  The  mathematical  models  of  the  batteries  are  based  on 
stochastic  approaches  or  empirical  equations  to  predict  the  run 
time,  capacity  and  efficiency  of  the  batteries  [1],  Furthermore,  the 
mathematical  models  do  not  have  direct  relation  between  the 
battery  model  parameters  and  the  I—V  (current-voltage)  charac¬ 
teristics  of  the  batteries:  hence  they  have  limited  use  in  circuit 
simulations. 
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The  electric  circuit  based  or  equivalent  circuit  models  (ECM)  are  the 
most  intuitive  for  use  in  circuit  simulations.  These  models  are  usually 
based  on  Thevenin  equivalents  and  impedances.  These  models  are 
simple  and  therefore  are  computationally  less  intensive.  One  of  the 
main  drawbacks  of  these  models  is  their  inability  to  predict  the  life¬ 
time  of  the  battery.  However,  their  simplicity  and  the  ability  to  predict 
the  I—V  characteristics  make  them  suitable  for  dynamic  modeling  of 
the  electric  vehicles.  An  ECM  comprises  of  a  combination  of  resistances 
and  capacitors.  These  parameters  are  multivariable  functions  of  SOC, 
current,  temperature,  number  of  charge  and  discharge  cycles  etc.  In 
order  to  use  ECMs  effectively  in  electric  vehicle  dynamic  simulation,  it 
is  critical  to  determine  the  parameters  of  the  circuit  accurately;  these 
parameters  are  usually  determined  by  extensive  experimentation  [5]. 
In  many  situations  the  evaluation  of  the  performance  of  an  electric 
drive  train  may  require  experiments  using  different  type  of  batteries 
from  different  manufacturers.  These  experiments  are  generally  very 
resource  intensive.  Therefore  it  is  generally  not  feasible  to  perform 
experiments  for  each  type  of  battery  to  extract  parameters  for  ECM.  To 
overcome  this  problem,  recently  Kumar  and  Bauer  [6]  have  suggested 
the  application  of  a  multi-objective  genetic  algorithm  (MOGA)  for 
extracting  the  battery  parameters. 

In  this  paper,  we  employ  a  modified  version  of  the  approach 
outlined  in  Ref.  [6]  for  extracting  the  battery  parameters  and 
determine  the  charge  and  discharge  characteristics  of  five  batteries 
(three  industry  manufactured  batteries  and  two  batteries  in 
research  and  development  stage).  The  genetic  algorithm  holds 
greater  benefits  compared  to  other  traditional  parameter- 
estimation  approaches  when  analytic  solutions  of  the  problem  do 
not  exist  and  when  the  number  of  design  parameters  are  large 
which  in  this  work  are  31.  The  genetic  algorithm  is  also  well-known 
for  its  minimal  requirements  of  prior  knowledge  of  the  solution¬ 
searching  optimization  problem.  The  computed  results  are 
compared  with  the  catalog  data  (in  case  of  industry  manufactured 
batteries)  and  the  experimental  data  (in  case  of  batteries  in 
research  and  development  stage). 

2.  Methodology 

Since  the  equivalent  circuit  model  is  described  by  a  set  of  non¬ 
linear  equations  [5,6],  it  is  not  tractable  for  obtaining  an  analyt¬ 
ical  solution.  However,  the  computations  employing  the  solution 
searching  technique  can  be  used  to  determine  the  parameter  values 
needed  for  modeling  the  batteries.  Genetic  algorithm  (GA)  appears 
to  be  an  ideal  candidate  for  such  task.  In  this  work,  a  modified 
version  of  multi-objective  non-dominated  sorting  genetic  algo¬ 
rithm  (NSGA-II)  has  been  employed  [7],  The  multi-objective  genetic 
algorithm  (MOGA)  code  was  first  validated  by  extracting  the  bat¬ 
tery  parameters  for  the  battery  [8]  used  in  Ref.  [6];  the  results  were 
compared  with  those  obtained  in  Ref.  [6],  After  the  validation,  the 
battery  modeling  code  was  employed  to  extract  the  battery  pa¬ 
rameters  of  two  additional  industrial  lithium-ion  batteries  [9,10] 
and  the  results  were  compared  with  the  catalog  data.  Finally  the 
code  was  applied  to  extract  the  optimal  parameters  of  recently 
proposed  high  performance  batteries,  three  Li— O2  batteries  [11] 
and  two  alkaline  batteries  with  an  iron  electrode  [12]  and  Potas¬ 
sium  Hydroxide  electrolyte;  these  two  batteries  are  in  early  stage  of 
research  and  development.  For  all  the  five  cases,  charging  and 
discharging  characteristics  were  obtained  which  were  compared 
with  either  the  catalog  data  (for  industrial  batteries)  or  with  the 
experimental  data  (for  research  domain  batteries). 

2.1.  Equivalent  circuit  battery  model 

The  schematic  of  the  equivalent  circuit  model  (ECM)  used  in  this 
work  is  shown  in  Fig.  1  [5,6],  The  parameters  in  ECM  are  identified 


Ci 


as  the  resistances  (Ri,  R2),  capacitance  (Ci),  and  open  circuit  voltage 
of  the  battery  (V0).  The  presence  of  one  pair  of  capacitance  and 
resistance,  i.e.  Ci  and  R2,  makes  the  model  to  be  linear  with  time. 
However,  non-linear  ECM  with  more  than  one  capacitance-resis¬ 
tance  pairs  has  also  been  proposed  [13,14],  which  is  supposed  to  be 
more  accurate  in  capturing  the  transient  behavior  of  the  battery. 
Non-linear  ECM  inevitably  introduces  more  modeling  parameters 
due  to  the  increased  model  complexity.  For  simplicity,  only  the 
linear  ECM  shown  in  Fig.  1  is  employed  in  this  paper. 

These  parameters  are  all  functions  of  the  current  State-of- 
Charge  (SOC)  and  discharge  rate  (Cr)  of  the  battery,  which  can  be 
formulated  as  following  polynomials  to  account  for  nonlinear 
phenomenon  in  the  battery  [6]: 

Resistance: 

R]  =  (a3  +  a2Cr  +  a3C2  )e~fl7xS0C  +  (a4  +  a5Cr  +  a6C2)  (1) 

R2  =  (as  +  a9Cr  +  aioCr2)e~a,4xSOC  +  (an  +  ai2Cr  +  a13C2) 

(2) 

Capacitance: 

Ci  =  -(ai5  +  a16Cr  +  a17C2je_a2lxSOC-i-  (ai8  +  a1gCr-i-a2oC2j 

(3) 

Open  circuit  voltage: 

Vq  =  —  (a22  +  a3oCr  +  a31  C2^e~a23xSOC  +  (a24  +  a25  x  SOC 
+  a26  x  SOC2  +  a27  x  SOC3  j  -  a28Cr  +  a29C2 

(4) 

In  Equations  (l)-(4),  the  coefficients  ai  through  a3i  are  model 
parameters  intrinsically  determined  by  the  battery  and  its  opera¬ 
tional  condition,  such  as  the  operation  temperature. 

The  time  dependent  voltage  output  of  the  battery  cell  is  then 
given  by  the  equation  [15]: 

V(t)  =  9^91e-tMQ  +  Vq  _  /Rj  -  ir2  (1  _  e-t/R2Cr)  (5) 

where  Qo  is  the  nominal  capacity  of  the  battery  cell  (A  h),  t  is  time 
duration  since  charge/discharge  starts  (seconds),  and  I  is  the  cur¬ 
rent  (A). 

Equations  (1  )-(5)  describe  the  performance  of  a  battery  by  V(t) 
as  a  function  of  SOC,  which  varies  from  0  to  100%  of  the  battery’s 
capacity.  Nominal  capacity  Qo  is  the  inherent  property  of  a  battery 
cell,  thus  it  remains  constant  and  can  be  measured.  Charging/ 
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Fig.  2.  Flowchart  of  the  optimization  process  for  determining  ai-a3i  using  the  modified  NSGA-II  MOGA. 


discharging  rate  Cr  is  controllable  based  on  the  demand  of  the 
application  and  is  likely  to  remain  constant  for  an  identical  oper¬ 
ation.  Therefore,  in  Equations  (1 )— (5)  the  model  parameters,  cq- 
03i,  remain  to  be  determined.  Once  the  values  of  01-031  are  known, 
the  complete  battery  model  is  established  and  its  charging/dis¬ 
charging  performance  can  be  predicted  by  Equations  (1)— (5).  This 
model  can  then  be  employed  in  simulating  the  performance  of  an 
electric  drive  train. 

3.  The  modified  multi-objective  NSGA-II  genetic  algorithm 

Computer  simulations  of  Equations  (1)— (5)  in  conjunction  with 
a  multi-objective  genetic  algorithm  (MOGA)  provide  an  effective 
approach  for  extracting  the  battery  parameters,  Oi— 031,  for  a 
particular  battery.  MOGA  is  used  to  find  the  values  of  these  co¬ 
efficients  cq—  <331  that  produce  a  voltage— SOC  characteristic  which 
best  fits  the  experimentally  measured  curve  by  the  battery 
manufacturer. 

Multi-objective  optimization  problems  require  the  optimization 
of  multiple  objectives.  However,  a  single  absolutely  optimal  solu¬ 
tion  satisfying  multiple  objectives  cannot  be  achieved  by  employ¬ 
ing  MOGA.  It  can  be  accomplished  by  a  genetic  algorithm  only  for  a 
single  optimization  objective.  Instead,  for  a  multi-objective  prob¬ 
lem  a  set  of  optimal  solutions  is  obtained  called  the  Pareto-front. 
The  area  which  the  Pareto-front  covers  is  called  the  hyper¬ 
volume  [16],  Any  point  inside  the  hyper-volume  is  considered  to 
be  dominated  by  at  least  one  solution  in  the  Pareto-front.  However, 
the  Pareto-front  itself  is  non-dominated;  meaning  no  one  solution 
in  this  set  can  optimally  satisfy  all  the  objectives. 

In  this  work,  non-dominated  sorting  genetic  algorithm  (NSGA- 
II)  is  used  [7],  The  code  based  on  NSGA-II  is  known  as  the  j-metal 
code  which  can  be  obtained  from  the  website  [17],  The  multi¬ 
objective  genetic  algorithm  (MOGA)  is  a  two-stage  process.  In 
this  process,  one  set  of  parameters  is  considered  to  be  one  indi¬ 
vidual.  First,  a  population  of  individuals  is  ranked  by  fitness  and  the 
best  are  selected  to  reproduce.  Reproduction  then  takes  place  with 
mutation  and  cross-over  to  create  an  intermediate  generation.  The 


new  generation  is  then  ranked  and  the  process  repeats  [18],  Mu¬ 
tation  randomly  occurs  at  low  probabilities  from  one  generation  to 
the  next,  altering  part  of  a  particular  individual.  Reproduction  is 
carried  out  through  the  random  pairing  of  two  individuals.  Once 
matched,  a  random  point  on  the  pair  is  chosen  where  cross-over 
occurs  and  the  two  strings  swap  chunks  of  code  [19].  The  selec¬ 
tion  process  is  carried  out  by  “binary  tournaments”.  A  population  is 
divided  into  groups  of  a  certain  size.  Each  group  competes  in  a 
tournament  in  which  the  individual  with  the  highest  fitness 
(lowest  objective  values)  wins  [20],  These  winners  are  chosen  as 
the  parents  for  the  next  generation  and  the  entire  process  is 
repeated. 

To  apply  MOGA  to  our  specific  problem  of  battery  parameter 
extraction,  standard  NSGA-II  code  [17]  was  modified.  Sixteen 
sample  points  were  read  from  the  catalog  voltage— SOC  charging/ 
discharging  characteristics  of  a  particular  battery.  Typically,  two 
exponential  regions  (at  high  and  low  SOC)  and  one  linear  region 
(between  high  and  low  SOC)  can  be  identified  for  most  voltage- 
SOC  characteristics.  Therefore,  it  is  best  to  have  more  sample  points 
in  the  exponential  regions  to  accurately  capture  the  steep  changes 
in  the  charging/discharging  characteristics,  and  keeping  fewer 
points  for  the  linear  region  in  order  to  save  the  computational 
resources. 


Table  1 

Batteries  considered  and  the  type  of  optimization  applied  to  extract 
parameters  Qi— 

Battery  Discharging 

e-PLB  [8]  M 

Lithium-ion  (ANR26650)  [9]  M 

Lithium-ion  (Seiko  MS614SE)  [10]  M 

U— 02  [11]  S 

Li-02  (+composite  carbon  cathode)  [11]  S 

Li-02  (+gold  loaded  carbon  cathode)  [11]  S 

Alkaline  battery  with  iron  electrode  [12]  S 

Alkaline  battery  with  iron  electrode  with  S 

bismuth  sulfide  ADDITIVE)  [12] 

*S  =  single  objective,  M  =  multi-objective. 
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Fig.  4.  Discharging  characteristics  for  Cr  =  0.5  C  (left)  and  4  C  (right)  for  ePLB  battery. 

where  Vc  is  the  catalog  voltage  value  and  Vg  is  the  model  output 
voltage  value  computed  from  the  genetic  algorithm  by  determining 
Qi—  03i  at  each  generation;  n  is  the  number  of  sample  points. 

Equation  (6)  describes  a  single-objective  optimization  problem. 
However,  recalling  that  the  voltage  characteristics  are  affected  by 
the  charging/discharging  rate,  a  multi-objective  optimization 
problem  is  formulated  as  follows: 

Objectives  1—  m  are  defined  by  the  functions: 

f  m  =  \Vc  -  ^kcr,  +  |Vc  -  Vg|2,Cri  +  ...  +  |Vc  -  Vg|niCri 

f2(x)  m  JVc  -  vg|1Cr2  +  \Vc  -  vg\2Cr2  + ...  +  \Vc  -  vg\nCri 

\fm(x)  m IVc  -  vg\XCrm  +  \vc  -  vs\2Crm  + |vc  -  vs\nCrm 

(7) 

where  m  is  the  number  of  different  charging/discharging  rates. 

The  objective  functions  in  Equation  (7)  are  minimized  by  using 
MOGA.  Closer  to  zero  these  functions  become,  better  agreement  is 
obtained  between  the  catalog  curve  and  the  model  output  curve. 
Since  negative  voltage  values  are  not  expected,  a  constraint  is 
implemented.  For  any  set  01—031  which  produces  negative  voltages, 
a  constraint  is  implemented  based  on  how  many  voltages  are  below 
zero  and  how  negative  the  values  of  these  voltages  are.  Individuals 
with  no  constraint  are  considered  better  fit  than  the  individuals 
with  a  constraint.  The  MOGA  simulation  is  considered  converged 
when  the  desired  minimum  for  the  objectives  in  Equation  (7)  is 
achieved.  The  converged  values  of  01—031  provide  the  closest  fit  to 
the  charging/discharging  curve  obtained  by  the  battery  model 


Since  the  ultimate  goal  is  to  determine  the  values  01-031  that 
would  lead  to  minimum  difference  between  the  model  output  and 
the  manufacturer  provided  battery  characteristic  curves,  the 
objective  function  is  defined  as  the  absolute  value  of  the  difference 
between  the  model  output  and  the  catalog  value  at  each  sample 
point.  The  objective  function  is  defined  as: 

Objective  :/(x)  =  \VC  -  Vg\r  +  |VC  -  Vg|2  +  ...  +  \VC  -  Vs\n 

(6) 


SOC  (%) 


Fig.  5.  Charging  characteristics  for  C,  =  2  C 
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Fig.  6.  Discharging  characteristics  for  Cr  =  0.33  C  (left)  and  10  C  (right). 


given  by  Equations  (1 )— (5),  when  compared  to  the  catalog  curve.  A 
flowchart  explaining  the  entire  optimization  process  for  deter¬ 
mining  ai-d3i  is  shown  in  Fig.  2. 

4.  Simulation  results 

In  this  paper,  using  MOGA,  parameters  are  extracted  for  eight 
batteries,  including  three  lithium-ion  batteries  (ePLB-cell  [8], 
MS614SE-cell  [9],  and  ANR26650Ml-cell  [10]),  three  Li-Oz  batte¬ 
ries  [11],  and  two  alkaline  batteries  with  iron  electrode  [12],  Based 


SOC(%) 


Fig.  9.  Charging  characteristics  for  5000  mA  g  1  current  density. 

on  the  data  availability  of  battery  characteristics  under  different 
charging/discharging  rates,  either  a  single-objective  optimization 
or  a  multi-objective  optimization  was  applied.  A  list  of  batteries 
considered  and  the  type  of  optimization  employed  is  given  in 
Table  1. 

4.1.  Battery  #1:  ePLB  cell  [8] 

The  ePLB-cell  [8]  was  modeled  using  ECM  and  MOGA  in  Ref.  [1  ]; 
this  case  has  been  used  for  our  code  validation  using  ECM  and 
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Fig.  8.  Discharging  characteristics  for  Cr  =  1  C  (left)  and  6.7  C  (right). 
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Fig.  11.  Charging  characteristics  for  70  mA  g-1  current  density. 


MOGA.  Two  charging  and  discharging  rates  of  Cr  =  0.5  C  and  4  C 
(C=  8  A)  were  considered.  The  extracted  battery  parameters  01-031 
are  summarized  in  the  appendices  in  Table  Al.  The  computed 
charging  and  discharging  characteristics  and  their  comparisons 
with  catalog  data  are  shown  in  Figs.  3  and  4  respectively. 

4.2.  Battery  #2:  ANR26650  cell  [9] 

ANR26650  battery  is  the  next  generation  of  A123  Systems  pio¬ 
neering  26650  cylindrical  cells.  It  is  a  high  power  battery  cell  with 
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Fig.  12.  Discharging  characteristics  for  70  mA  g  1  current  density. 


large  energy  density  and  lower  impedance.  Due  to  its  versatility, 
the  ANR26650  battery  is  used  in  a  wide  variety  of  applications  e.g. 
in  EVs  and  HEVs.  Fig.  5  shows  the  comparison  of  our  model’s  output 
with  catalog  data  [9]  for  charging  characteristics  for  Cr  =  2  C 
(C  =  3  A)  and  Fig.  6  shows  the  comparison  of  our  model’s  output 
with  catalog  data  [9]  for  discharging  characteristics  for  Cr  of  0.33  C 
and  10  C  (C  =  3  A).  Table  Al  in  appendices  summarizes  the 
extracted  parameters  01—031  for  this  battery. 

4.3.  Battery  #3:  MS614SE  cell  [10] 

MS614SE  cell  is  a  high-performance  rechargeable  micro  battery 
produced  by  Seiko  Instruments  for  sophisticated  electronic  in¬ 
struments  and  equipment.  Fig.  7  shows  the  comparison  of  our 
model’s  output  with  catalog  data  [10]  for  charging  characteristics 
for  Cr  =  6.7  C  (C  =  15  pA)  and  Fig.  8  shows  the  comparison  of  our 
model’s  output  with  catalog  data  [10]  for  discharging  characteris¬ 
tics  for  Cr  of  1  C  and  6.7  C  (C  =  15  pA).  Table  Al  in  appendices 
summarizes  the  extracted  parameters  01—031  for  this  battery. 

4.4.  Battery  #4:  Li-02  cell  [11] 

In  a  recent  research  paper,  rechargeable  non-aqueous  Li— O2 
battery  has  been  proposed  [11];  it  is  receiving  a  great  deal  of  in¬ 
terest  because  of  its  high  specific  energy  which  far  exceeds  the  best 
achievable  by  lithium-ion  cells.  Fig.  9  shows  the  comparison  of  our 
model’s  output  with  the  experimental  data  [11  ]  for  charging  char¬ 
acteristics  under  5000  mA  g-1  charging  current  density 
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Fig.  14.  Discharging  characteristics  for  70  mA  g-1  current  density. 


!  et  at  /  j 


735 


SOC(%)  SOC(%) 


100  80  60  40  20  0 

SOC  (%) 

Fig.  16.  Discharging  characteristics  for  Cr  =  0.05  C. 

(normalized  by  the  mass  of  carbon)  and  Fig.  10  shows  the  com¬ 
parison  of  our  model’s  output  with  data  [11]  for  discharging  char¬ 
acteristics  under  5000  mA  g  1  discharging  current  density 
(normalized  by  the  mass  of  carbon).  The  extracted  battery  param¬ 
eters  01—031  are  summarized  in  the  appendices  in  Table  Al. 

4.5.  Battery  #5:  U—O2  cell  (with  composite  carbon  cathode)  [11] 

Figs.  11  and  12  show  the  comparison  of  our  model’s  output  with 
the  experimental  data  [11]  for  charging  and  discharging  charac¬ 
teristics  under  70  mA  g-1  current  density  (normalized  by  the  mass 
of  carbon)  respectively.  The  extracted  battery  parameters  Oi— 031 
are  summarized  in  the  appendices  in  Table  Al. 

4.6.  Battery  #6:  U—O2  cell  (with  gold  loaded  composite  carbon 
cathode)  [11] 

Figs.  13  and  14  show  the  comparison  of  our  model’s  output  with 
the  experimental  data  [11]  for  charging  and  discharging  charac¬ 
teristics  under  70  mA  g-1  current  density  (normalized  by  the  mass 
of  carbon)  respectively.  The  extracted  battery  parameters  01-031 
are  summarized  in  the  appendices  in  Table  Al. 

4.7.  Battery  #7:  alkaline  battery  with  iron  electrode  and  KOH 
electrolyte  cell  [12] 


and  large-scale  implementation  of  energy  storage.  Batteries  with 
carbonyl  iron  electrodes  have  been  studied  in  Ref.  [12]  where  their 
charging  and  discharging  characteristics  have  been  reported.  Fig.  15 
shows  the  comparison  of  our  model’s  output  with  the  experimental 
data  [12]  for  charging  characteristics  for  Cr  =  0.5  C  (C  =  0.36  A)  and 
Fig.  16  shows  the  comparison  of  our  model’s  output  with  data  [12] 
for  discharging  characteristics  for  Cr  =  0.05  C  (C  =  0.36  A).  The 
extracted  battery  parameters  01—031  are  summarized  in  Table  Al. 

4.8.  Battery  #8:  alkaline  battery  with  iron  electrode  with  bismuth 
sulfide  and  KOH  electrolyte  cell  [12] 

To  improve  the  specific  capacity  and  charging  efficiency  of  iron 
batteries,  bismuth  sulfide  can  be  added  [12],  Figs.  17  and  18  show 
the  comparison  of  our  model’s  output  with  the  experimental  data 
[12]  for  charging  characteristics  for  Cr  =  0.5  C  (C  =  0.64  A)  and 
discharging  characteristics  for  Cr  =  0.05  C  (C  —  0.64  A)  respectively. 
The  extracted  battery  parameters  ai— 031  are  summarized  in 
Table  Al. 


5.  Conclusions 


In  this  work,  a  modified  NSGA-I1  multi-objective  genetic  algo¬ 
rithm  [7]  based  optimization  code  has  been  developed  to  effec¬ 
tively  perform  the  extraction  of  battery  parameters  required  for 
modeling  the  charging/discharging  performance  of  a  battery  using 
an  equivalent  circuit  model  for  application  in  electric  drive  train 
simulation  of  electric  vehicles.  The  results  obtained  from  our  code 
for  eight  different  batteries  are  in  reasonably  good  agreement  with 
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Recently,  iron-based  alkaline  batteries  have  been  identified  as  a 
potential  solution  in  addressing  the  challenges  of  durability,  cost, 


Fig.  18.  Discharging  characteristics  for  Cr  =  0.05  C. 
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the  catalog/experimental  data  for  charging/discharging  character¬ 
istics.  The  results  suggest  that  the  developed  code  can  be  used  with 
confidence  for  extraction  of  battery  parameters  needed  in  an  ECM 
model  for  a  large  variety  of  batteries  with  different  sizes,  capacities, 
materials,  etc.  It  should  however  be  noted  that  the  model  outputs  of 


charging/discharging  characteristics  from  our  code  may  slightly 
differ  from  the  catalog  characteristics  in  some  cases,  especially  at 
low  and  high  SOC  when  drastic  changes  in  voltage  may  occur.  This 
suggests  future  improvement  in  the  code  may  include  weighting  of 
the  point  distribution  in  these  regions. 


Appendices 


Table  A1 

Battery  parameters  di—  CZ31  for  various  batteries  under  consideration 


Validation  ANR26650 

Discharging  Charging  Discharging 


°21 


025 

026 

027 

028 


0.20835291 

0.03579959 

0.00039606 

0.15919108 

0.0058729 

0.014583 

0.84904339 

0.01799999 

0.07685504 

0.4562723 

3.9482E-05 

6.3776E-07 

2.2766E-06 

38.8827613 

0.01402422 

0.13564355 

0.01853951 

9.92993734 

9.708302 

8.79748508 

0.67106684 

1.69655185 

0.52269596 

1.69347882 

0.27104638 

0.13853725 

0.33570299 

0.86261017 

0.73759611 

0.81575463 

0.71762382 


0.30689735 

0.01562162 

0.00685809 

0.10903136 

0.01125884 

0.0124507 

0.13708388 

5.51  E-05 

1.38E-05 

5.10E-06 

1.92E-05 

1.68E-05 

3.47E-06 

0.92172096 

0.45025886 

0.37415859 

0.03503156 

0.90074554 

0.83842829 

0.9166211 

0.5192527 

0.03956599 

0.12123495 

3.57513804 

0.01323247 

3.73E-05 

0.17322835 

0.61028401 

0.72975429 

0.34930874 

0.9163944 


0.00461877 

4.8659E-06 

0.01179723 

0.7512944 

0.07957936 

0.017111 

0.37531505 

0.89502179 

0.45275245 

0.93659989 

0.00071951 

0.0013407 

0.03375708 

45.9334064 

1.03289541 

9.67102855 

7.67531457 

9.96331567 

6.62020206 

9.91798063 

3.39887702 

2.18565757 

0.297832 

1.42529313 

0.34651381 

0.01210208 

0.54727467 

0.17518266 

0.6578156 

0.75252781 

0.47158002 


Li-02  carbon 


Li-02  gold  loaded 


Seiko  MS614SE 

Discharging 

0.091205243 

0.540474998 

0.820580154 

0.118955849 

0.946536 

0.803583273 

0.773938703 

8.628790134 

1.691575926 

0.674445879 

0.036847234 

0.994080133 

0.918395111 

17.14106692 

-70.2324694 

-68.40602547 

7.37606324 

16.54940513 

1.424637761 

21.5736066 

7.91293592 

1.249609558 

0.160405955 


0.000793251 

1.91075E-05 

5.272873439 

0.876496134 

0.10467427 

0.873173876 


Li— 02 


Discharging 

0.01073646 

0.18975732 

0.78228093 

0.17712076 

0.97380714 

0.57579603 

0.96127054 

4.57E-05 

0.29136878 

0.01467387 

8.99E-06 

1.51  E-05 

3.14E-06 

12.2351843 

0.43140728 

0.53229058 

0.32551931 

9.18060086 

9.92182253 

9.76708248 

0.23486365 

0.50907584 

0.60452313 

1.97255537 

0.66879075 

8.34E-04 

0.02384203 

0.15386122 

0.34071846 

0.80860882 

0.34129681 


Charging 

53.5953143 

1.42771541 

10.5988734 

0.05123828 

0.16838811 

0.53580137 

19.7204805 

0.00460395 

0.62259064 

0.02035467 

0.14779553 

0.79662184 

0.95324255 

69.9632914 

0.31602659 

0.59819497 

0.14692834 

4.14873414 

2.92392616 

7.15771813 

0.29390599 

2.38881108 

0.19900778 

1.90399394 

2.13001815 

4.37E-04 

0.00190092 

0.78542756 

0.11985901 

0.80790961 

0.38748586 


Iron  +  bismuth  sulfide 


Discharging  Charging  Discharging  Charging  Discharging  Charging  Discharging  Charging 


a\ 

a2 


0.82432906 

0.16163508 

0.99902869 

0.79862881 

0.64920124 

0.41082605 

0.53770096 

0.99676323 

0.98328136 

0.86287458 

1.93E-05 

4.35E-04 

0.00321498 

15.019207 

0.07644745 

0.33011917 

0.28219357 

9.88748962 

9.99319013 

8.42971362 

0.48624695 

1.29634101 

0.20261211 

1.46582257 

0.00295324 

0.19549466 


14.5241027 

24.4756143 

2.63245697 

0.25224944 

0.1630204 

0.00278872 

119.108559 

0.99948113 

0.99956374 

0.93117929 

0.02263178 

0.4718394 

0.08762695 

13.3261146 

0.09699596 

0.94621908 

0.31727155 

0.00364755 

0.82662482 

7.16279946 

0.17460783 

1.20141209 

1.40E-04 

2.57245923 

0.5439663 

1.69E-04 


0.26549751 

0.66016791 

0.93141066 

0.27658709 

0.15746281 

0.49241206 

0.48588926 

0.61531353 

0.46278238 

0.21019478 

5.46E-06 

5.23E-05 

0.00705363 

14.8134438 

0.04481933 

0.29095625 

0.62979428 

9.54738307 

9.05491208 

9.88223081 

0.93659762 

1.5572842 

0.10064423 

1.00973436 

0.09635125 

0.11227219 


14.1224938 

4.40953222 

37.0450311 

0.00590928 

0.3036228 

0.56180101 

75.4897042 

0.00317614 

0.28678398 

0.0982417 

0.0382265 

0.31200386 

0.30672593 

144.115905 

0.60145006 

0.44992672 

0.50843293 

0.19036418 

3.52250087 

9.60579147 

0.06051375 

2.76280474 

0.66444894 

0.85216059 

2.10345909 

0.33182808 


29.4914986 

72.9081168 

73.0177209 

0.26222204 

0.09864263 

0.78515787 

149.996207 

0.08526336 

0.25717061 

0.646143 

5.60E-05 

3.20E-04 

0.02603126 

149.983968 

3.41  E-06 

0.04565322 

0.09046255 

9.99996127 

9.97043314 

7.92405725 

0.85690665 

0.63327953 

0.23745515 

0.29576144 

0.19194707 

6.31E-07 


0.45462307 

0.54913924 

0.12070674 

0.75297676 

0.98736132 

0.11234894 

0.30099129 

8.61  E-ll 

1.09E-07 

4.69E-06 

3.33E-07 

3.08E-08 

4.01  E-10 

0.4571012 

0.06429345 

0.25185158 

0.39543831 

0.92995916 

0.83240631 

0.9479404 

0.8306388 

0.08588711 

0.26704514 

2.37981823 

0.009514 

0.1029671 


14.5599838 

11.1430598 

38.1558624 

0.25523467 

0.81490269 

0.61007881 

88.931616 

0.85858184 

0.06270029 

0.87188255 

8.11E-08 

0.00160551 

0.02568106 

149.852995 

0.02258402 

0.01134694 

0.29202957 

9.90334007 

8.35353122 

7.82615417 

0.20653592 

0.02121622 

0.14166353 

0.93034068 

0.01259381 

0.01995769 


0.55333126 

0.9640732 

0.40420978 

0.54544581 

1.95806029 

0.73571603 

0.17610559 

1.08E-10 

2.76E-10 

6.66E-08 

1.12E-09 

6.40E-09 

2.39E-09 

0.06898583 

0.00345412 

0.00376443 

0.23512432 

8.71964227 

7.96101789 

8.82767317 

0.65028507 

0.50355483 

0.16544309 

0.39883709 

0.10754302 

0.00102557 
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Table  Al  ( continued ) 


Li— 02  carbon  Li-02  gold  loaded  Iron  Iron  +  bismuth  sulfide 

Discharging  Charging  Discharging  Charging  Discharging  Charging  Discharging  Charging 


0.02204306 

0.29079774 

0.88948186 

0.26938344 

0.05713941 


3.47E-04 

0.81658536 

0.84944614 

0.64625688 

0.35704537 


3.84E-04 

0.83635088 

0.04848865 

0.91331299 

0.83708787 


0.5021228 

0.82468851 

0.507183 

0.05154254 

0.00115925 


2.74E-07 

0.9616469 

0.4414678 

0.48510127 

0.44912057 


0.00852127 

0.35565489 

0.0153644 

0.52007456 

0.25172147 


0.0052646 

0.54176887 

0.16401757 

0.36083715 

0.19398491 


0.14710023 

0.27629373 

0.81392545 

0.65849334 

0.96059927 
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